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Abstract: We study the effects of strict conservation laws and the problem of negative contributions to final 
momentum distribution during the freeze out through 3-dimcnsional hypersurfaces with space-like normal. We study 
some suggested solutions for this problem, and demonstrate it on one example. 
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■ 1. Introduction 

Fluid dynamical models, especially their simpler versions are very popular in heavy ion physics, because they con- 
nect directly collective macroscopic matter properties, like the Equation of State (EoS) or transport properties, to 

■ measurables. 

. ^ | Particles which leave the system and reach the detectors, can be taken into account via source (drain) terms in 
the 4-dimensional space-time based on kinetic considerations, or in a more simplified way via freeze out (FO) or final 
break-up schemes, where the frozen out particles are formed on a 3-dimensional hypersurface in space-time. This 
information is then used as input to compute measurables such as two-particle correlation, transverse-, longitudinal-, 
radial-, and cylindrical- flow, transverse momentum and transverse mass spectra, etc. 

In this paper we concentrate on freeze out. A basic standard assumption in this case is that freeze out happens 
across a hypersurface as already mentioned, so it can be pictured as a discontinuity where the kinetic properties of 
the matter, such as energy density and momentum distribution change suddenly. The hypersurface is an idealization 
of a layer of finite thickness (of the order of a mean free path or collision time) where the frozen-out particles are 
formed and the interactions in the matter bacome negligible. The dynamics af-this layer is described in different 
kinetic models such as Monte Carlo modelsErB or four- volume emission modelsa Tl In fact, the zero thickness limit 
of such a layer is an over-idealization of kinetic freeze out in heavy ion reactions, while it is applicable on more 
macroscopic scales like in astrophysics 0. 



*On the other hand, if kinetic freeze out coincides with a rapid phase transition, like in the case of rapid deconfinement 
transition of supercooled quark-gluon plasma, the sharp freeze out hypersurface idealization may still be applicable even for 
heavy ion reactions. It is, however, beyond the scope of this work to study the freeze out dynamics and kinetics in this latter 
case. 
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Two types of hypersurfaces are distinguished: those with a space-like normal vector, da^da^ = —(da) 2 (e.g. 
events happening on a propagating 2-dimensional surface) and those with a time-like normal vector da^da „ = (da) 2 
(a common example of which is an overall sudden change in a finite volume) . 

Once the freeze out surface is determined, one can compute measurables. Landau when drafting his hydrodynamical 
modefl, just evaluated the flow velocity distribution at freeze out, and this distribution served as a basis, for all 
observables. This approach was used in early fluid dynamical simulations of heavy ion collisions, .alsomli This 
procedure was improved to add thermal velocities to the flow velocities at freeze out, by-MilekhinExEj and later by 
Cooper and FryeuJ. This method is widely used, however it rises at least three problemsE-3. 

First, in some cases before the 90's, the possible existence of discpatiHuities across hypersurfaces with time-like 
normal vectors was not taken into account or considered unphysicaFpn23. This point was studied recentlya so we 
do not discuss it further. 

Second, since the kinetic properties of the matter are different on the two sides of the front, the explicit evaluation 
of conservation laws aecoss^the freeze out surface should be taken into account which is not alwaws easy to implement. 
In some (simple) casesBjLJ, these conservation laws are enforced and discussed. For example inc3, it was pointed out 
that the freeze out momentum distribution for hypersurfaces with time-like normal may become locally anisotropic. 
InnZl a solution for post FO massless, baryonfree Bose gas was presented. Here we remind the procedure that should 
be followed in section ||. 

The third problem is a conceptual problem arising in the Cooper- Frye freeze out description when we apply it to 
a hypersurface with space- like normal: it is the problem of negative contributions (see section^). This is the main 
subject of this paper. This problem appears in all freeze out calculations up to now we are aware of, and to our 
knowledge it was not satisfactorily discussed .yet .in. the literature. It was recognized by some of those who apalied 
the Cooper-Frye freeze out description beforec3c3'C3. A possible partial solution was presented ia part 2 of refHj for 
noninteracting massless particles, in 1+1 dimension using the post FO cut Jiittner ansatz. Inc3 it was shown that 
in an oversimplified kinetic freeze out model one can .obtain the cut Jiittner distribution as post FO distribution. 
In section [| we complete and generalize the results oin3 and present an example for the solution of the freeze out 
problem. In section |g we suggest improvements that go beyond the cut Jiittner ansatz. 



2. Conservation laws 

across idealized freeze out discontinuities 

In the zero width limit of the freeze out domain (freeze out surface), the energy - momentum tenso ij | Chan ges dis- 
continuously across this surface. Consequently, the four-vector of the flow velocity may also change£ZrEjl23. These 
changes should be discussed in terms of the conservation laws. 

The invariant number of conserved particles (world lines) crossing a surface element, da 11 , is 

dN = da^ , (1) 

and the total number of all the particles crossing the FO hyper-surface, S, is 

N = ( N» da^ . (2) 
Js 

If we insert the kinetic definition of 

N» = J^-P 11 f FO (x,p) , (3) 



discussed discontinuities across propagating hypersurfaces, which have a space-like normal vector. If one applies 
Taub's formalism from 1948, to freeze out surfaces with time-like normal vectors, one gets a usual Taub adiabat but the 
equation of the Rayleigh line yields imaginary values for the particle current across the front. Thus, these-Jnypersurfaees were 
thought unphysical. However more recently, Taub's approach has been generalized to these hypersurfacest-J (see alscu3) while 
eliminating the imaginary particle currents arising from the equation of the Rayleigh line. Thus, it is possible to take into 
account conservation laws exactly across any surface of discontinuity with relativistic flow. 
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a0: 



into eq. ([I]) we obtain the Cooper-Frye formul. 

E W~ = / Ifo{x ' p) P ^ ' (4) 

where Jfo(XtP) is the post FO phase space distribution of frozen-out particles which is not known from the fluid 
dynamical model. The problem is to choose its form correctly. Usually one assumes that the pre-FO momentum 
distribution as well as the post FO distribution are both local thermal equilibrium distributions, with the same 
temperature, boosted by the local collective flow velocity on the actual side of the freeze out surface, although 
the post FO distribution need not be a thermal distribution. Parametrizing the post FO distribution as thermal, 
fFoiXjP'jTjn^u^), where T does not necessarily coincide with the pre FO temperature, and knowing the pre-FO 
baryon current and energy.- momentum tensor, Nq and Tq , we can calculate the post freeze out quantities N 11 and 
T^ v from the relationMrJ 

[N" = and [T» v da^} = 0, (5) 

across a surface element^] of normal vector da' 1 . Here [A] = A — Aq. The post FO distribution is not a thermal 
equilibrium distribution, so temperature does not exist, nevertheless, the conservation laws fix the parameters, e.g. 
T, n, u M , of our momentum distribution, fpo{x 1 p\T,n,u u ). 

To obtain a physically realizable result, in addition we have to check the condition for entropy increase: 

[S» or fl = |^i>l, (6) 

where, for both equilibrium and nonequilibrium FO distributionslll 
d 3 p 



p" f F o(x,p) [ln{(2^) 3 f FO (x,p)} - l] . (7) 



This condition is not necessary to obtain a solution of the freeze out problem, but it should always be checked 
to exclude non-physical solutions. We have to note at this stage, that the post FO distribution must not be an 
equilibrium (or stationary) solution of the Boltzmann Transport Equation, and consequently on the post FO side 
the energy flow-, baryon flow- and entropy flow- velocities may all be different. 

We can now remind briefly what the problem of negative contributions to the Cooper-Frye formula is, and a 
possible way out. For a FO surface with time-like normal, both p M and g?ct m are time-like vectors, thus, 

pPdap > 0, 

and the integrand in the integral (Q) is always positive. For a FO surface with space-like normal, p M is time-like and 
da^ is space-like, thus, p^da^ can be both positive and negative. (Note that p^ may point now both in the post- 
and pre- FO directions.) Thus, the integrand in the integral (|J) may change sign in the integration domain, and this 
indicates that part of the distribution contributes to a current going back, into the front while another part is coming 
out of the front. On the pre-FO side p M is unrestricted and p^da^ may have both signs, because we are supposing 
that pre-FO phase is in the thermal equilibrium. However, in the zero width limit of the FO front, it is difficult to 
understand such a situation. What happens actually is that internal rescatterings occur inside the finite FO domain 
and feed particles back to the pre-FO side to maintain the thermal equilibrium there. On the post-FO side, however, 
we do not allow rescattering and back scattering any more. If a particle has, passed the freeze out domain it cannot 
scatter back. In other words, the post-FO distribution should have the forir£3'E3, 

f FO (x,p,da») = f FO {x,p)Q{p» da^) , (8) 



*In numerical calculations the local freeze out surface can be determined most accurately via self-consistent iteratio 
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where 0(x) is the step function. Consequently, this distribution cannot be an ideal gas distribution. (On the pre-FO 
side, the distribution may or not be ideal). The conservation laws across a small element of the freeze out front with 
space-like normal take the form: 



ai 



d \ f* F o(^P,d^)p^) da, = I Ng(x) da, , (9) 



P 



^ f* FO (x,p,^)^ P v ) da,= l T^{x)da,. (10) 



P 



3. The allowed momentum region for space-like FO 

Let us assume that the FO process happens in the positive x direction, in other words we go in the positive x-direction 
from the pre FO domain to the post FO domain. The FO hyper-surface has a space-like normal (da^da, = —(da) 2 ), 
so that da, is orthogonal to the hyper-surface (i.e. to the time-like tangent vector of the surface i M , t^da, = 0) and 
points into the post FO (positive x- ) direction (while da^ points in the pre FO direction). 

Depending on the reference frame, the space-like FO front can propagate both in the positive or negative in- 
direction, or it can be Lorentz transformed into its own rest frame, the Rest Frame of the Front (RFF), where da, 
— (0, 1, 0, 0)da. In other reference frames da, = j(v, 1, 0, 0)da, where 7 = 1/ yl — v 2 . The parameter v = da$ /da x 
= —t x /t° is frame dependent and may be both positive and negative. 

In order to use the parameter v in the following discussion, we have to select a given reference frame and fix the 
value in that tame. Let us choose the frame comoving with the peak of the post FO invariant momentum distribution. 
Following ref.Ej let us denote this frame as the Rest Frame of the Gas (RFG). Note, however, that, contrary to what 
its name seems to suggest, this frame is not the Local Rest frame of the post FO matter, since the post FO distribution 
is not spherically symmetric in momentum space! Thus, for the following discussion we define v in the RFG frame 
as 

da 

v = T~ 

aa x RFG 

or we can also define and have the same w-value, by using the velocity of the peak of the post FO distribution, u RFG 
in the RFF. (Note: this is not equal with the post FO flow velocities, neither with the Eckart-, nor with the Landau- 
flow velocity!) Thus, 

„, _ U R.FG 





U RFG RFF 

and this means that in the RFF the peak four-velocity has the form u RFG = 7(1, v, 0, 0)\rff ■ This velocity v can 
also be both positive and negative (!) in RFF, i.e. the peak velocity may point to the post FO direction (as we would 
expect), and also in the pre FO direction in RFF which seems to be confusing. This is, however, not a problem in 
itself because both flow velocities (Eckart and Landau) are always positive on the post FO side. As we will see later, 
in special cases it is possible that we obtain negative post FO peak velocity for positive pre FO flow velocities. This 
indicates we have to discuss the importance of the cut by Q^p^da,), otherwise one might be tempted to believe that 
this cut is not affecting the post FO distribution too much, and the correct treatment causes only a few percent cut 
which is negligible. We show in the following that this is not the case. 

The p^da, > requirement in the RFG frame means that only momenta with component 



p x > -v^m 2 + (p*) 2 + {p±_) 2 , (11) 

contribute to the post FO momentum distribution. This means that the boundaries of the allowed domain in the 
[p x ,p±]- plane arc hyperbolas in the post FO RFG 
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f^j-(p^T=m\ (12) 

and the domains of the positive x- side of the corresponding hypcrboloids in the 3-dimcnsional momentum space 
may contribute to the FO distribution (Fig. [I]). In the massless limit the hyperboloids become cones around the 
x-axis and centered at the origin. 
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Fig. 1. The boundaries of the cut post FO distributions in the phase space indicating the regions where particles are allowed to 
freeze out. These boundaries are hyperboloids, for RFF velocities of: v = 0.65, v = and v = —0.35 (from left to right). Particles to the 
right of these hyperboloids may freeze out. While for large positive v values the cut is a small perturbation, for moderate or particularly 
for negative values of v the cut is far from negligible. The dotted lines show the massless limits. 



a. Dominant case. For v = the hyperboloid becomes a plane boundary at p x > in RFG, i.e. the boundary 
cuts the Jiittner distribution in the middle (because in the RFG the peak of the distribution is centered at p — 0). 
If v > 0, the peak velocity points to the post FO direction and the cut-hyperboloid is fully in the negative p x region, 
thus for large v values only a smaller fraction is excluded from the full Jiittner distribution. Most of the particles 
described by a full Jiittner distribution freeze out in this case. We will consider this as the dominant FO case. 

On the other hand, as we can see in Fig. |l|, for the v < case the peak flow points back into the pre FO direction, 
and only a small fraction of a full Jiittner distribution will freeze out. The cut eliminates the major part of the 
particles. This situation leads to less frozen out particles, but it yields an unusual post FO distribution. 



4. Conserved currents for cut Jiittner distribution 



We now study the particular case where Jfo is a Jiittner (or relativistic BoltzmannElH) distribution and so fp Q is 
a cut Jiittner distribution: 



f* FO (p) - f Juttn -(p)e^da ll ) = exp ( M P f RFG ) , (13) 
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where fi is the chemical potential related to the invariant scalar density n, of the non-cut Juttner distribution as fj, = 
Tin [n(27r?i) 3 / (4irm 2 T K 2 (y))] > an( i U RFG' " an( i ^ are parameters of the distribution f FO originating from the full 
Juttner distribution. These arc not the flow velocity, proper density and temperature of the cut Juttner distribution! 
The cut Juttner distribution is not a thermal equilibrium distribution, n^£. it does-jaot have a temperature at all. 

This distribution for massless particles was considered in part 2 of ref rJ (also irEj) . The cut selects particles with 
momenta p^da^ > 0. 

We can evaluate the baryon four-current, N^, by inserting the cut Juttner distribution into the definition, eq. 
(||), and we get a time directed-, N°, as well as a spatial component, N x , (where x is the direction of the spatial 
component of the FO-normal, da 11 , in RFG). In ref £3 the spatial component, N x , of the four current was not 
evaluated, so seeing only the zeroth component, N°, the unsuspecting reader might have believed falsly, that RFG 
(Rest Frame of the Gas) is the Local Rest frame of the gas. Performing the calculation, in the post FO RFG frame 
the baryon current reads as 

N° = f {vA + a 2 j [(l+j)K 2 (a) - IC 2 (a, b)} + 

j^V"} ^5 n(»,T)^, (14) 

N x = § [(l-v 2 )A-a 2 e- b ] n(p,T)^ , 

where j = sign(u), n — BttT 3, e^/ T (2ith)~ 3 , a — so that h(fj,, T) = na 2 K 2 {a)/2 is the invariant scalar density of 
the symmetric Juttner gas, b = a/yl — v 2 , v — dao/dcr x , A = (2 + 2b + b 2 )e~ b , and 

K n ( Z ,w) = ±g± z ~n dx (a .2_ z 2 )n -l/ 2 e -z 

i.e. K. n {z, z) = K n (z). Just as in case of the non-cut distributions the cut Juttner distribution yields few modified 
Bessel functions in the expression of the four currents, while the relativistic Fermi and relativistic Bose distributions 

lead to a series of these functions. When evaluating the limits we used the relation K n {a,b) ^> K n (a) 
2 n ~ 1 (n — l)!a _ ™. This baryon current may then be Lorentz transformed into the Eckart Local Rest (ELR) frame 
of the post FO matter, which moves with = /{N v N^) 1 / 2 = 7e(1, ve, 0, 0)\rfg in the RFG, or alternatively 
into the Rest Frame of the Freeze out front (RFF), where da^ = (0, 1, 0, 0)(da) 2 \fjFF and the velocity of the RFG is 
u E RFG = 7 ff (l, v , 0, 0)\rff ■ Then the Eckart flow velocity of the matter represented by the cut Juttner distribution 
viewed from the RFF is u E — 7 C (1, v c , 0, 0)\rff, where v c = (v + «e)/(1 + we)- 
The proper density (i.e. the density in the ELR frame) is obtained as 



n(fM,T,v) = ^/WW V . (15) 

Note that the proper density of the cut Juttner distribution, n, is reduced compared to the proper density of the 
complete spherical Juttner distribution, h. 

The energy momentum tensor in the post FO RFG is 



T oo = 3|T fiol [x 2 (a)+fif 1 (a)]-[/C 2 (a ) 6)+f/C 1 (a,6)]}+Bv}, 



rpOx _ 3hT 
1 4 



rpXX 



(l-v 2 )B-^-(b+l)e- b }, 
(jT [(WW - ^(a, b)} + v 3 B} , 



(16) 



3nT 
4 



[v(l-£)B + ^ [(l+j)K 2 (a) - tC 2 (a, b)} - ^f(b+l) e - b } , 

where B = (1 + b + b 2 /2 + b 3 /6)e~ b and T zz = T vv . This energy-momentum tensor may then be Lorentz transformed 
into the Landau Local Rest (LLR) frame of the post FO matter, which moves with u L in the RFG, or into the Rest 
Frame of the Freeze out Front (RFF) where da M = (0, 1, 0, Q)da. Alternatively both can be transformed to the frame 
where we want to evaluate the conservation laws, eq. (j|), and the parameters of the post FO, cut Juttner distribution 
can be determined so, that it satisfies the conservation laws. In the massless limit the energy momentum tensor in 
the RFG is: 

T 00 = 3nT (v+X) /2 , T 0x = 3nT (l-w 2 ) /4 , T xx — nT (v 3 +l) /2 



G 



and T zz = = (T 00 - T xx )/2. 

In addition we have to check the entropy condition. In the RFG frame the entropy current reads 

S° = | {(1-£M + QvB + (1 - T) fl2 J [(!+i)^2(a) - IC 2 (a, b)} + 

ja^l+m^-JC^b)}} ^ *i$H( w + i)( 4 _£) , 
S»= |[(l- v 2 )(l-^)A + 6(l-« 2 )B- a 2 (2 + 6-^)e- b ] ^ 

=S ^(l-^ 2 )(4-^) . 

Note that in the m = limit the vectors .S'' 1 and 7V M are parallel to each other. This is explained by Figure [j], which 
shows that in the RFG the cut in the m = limit becomes a central cone, and since the distribution is centrally 
symmetric in this frame, the integrals will be proportional to each other. 



A. Solubility of the freeze out problem 

The situation is non-trivial and we have to take into account the possible directions of the flow and of da^. Note: 
we must not assume that the flow is parallel to the freeze out direction. 

Let us start on the pre FO side labeled by '0'. Here in the LR frame Uq = (1, 0, 0, 0)\ pr e-LR and we can choose 
the x-direction in this frame to point into the FO direction, so that da^ — 70(^0, 1, 0, Q)da\ pre -LR- We assume that 
we know the FO hypersurface, i.e. we know vq. Then, in this frame the conservation laws have three nonvanishing 
components yielding three known parameters N^da^, T^da^, and T^da^. 

To find the solution we need these values in the RFG frame. However, the 3-dimensional direction of the x- 
axis will not change because the front is assumed to be isotropic in its own [y, z]-plane. Thus, in the RFG the 
peak flow parameter is u RFG — (1,0,0,0)1^0, an d the normal of the FO front is — j(v, 1,0, 0). Note that 
v =/= vq\ Furthermore, let us recall that the parameter v determines the post FO peak flow parameter in RFF, 
u rfg = v,0,0)\rff (where da^ = (M, 0, 0)da\ RFF ). 

Consequently the conservation laws (P,|lO|) yield three nonvanishing equations in the RFG frame, 

[iVMcrJ = , [T 0fi da^] = and [T^da^ = , 

which can determine the three unknown parameters of the post FO cut-Jiittner distribution, v, T, and n (or /i). 
While the first equation is an invariant scalar, the remaining two are components of a 4-vector, so they should be 
transformed into the same reference frame, i.e., to RFG. Since we evaluated the quantities based on the cut-Jiittner 
distribution in the RFG, we also need the pre FO quantities in the RFG. These can be determined by using the 
standard fluid dynamical form of T^ v as seen from the RFG. From this frame the pre FO flow velocity is given by 
the difference of the pre- and post-FO flow velocities: Uq = 7or(1, vqu, 0, 0)\rfg, where vqh = (vq —v)/(1 — vq v). 
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Fig. 2. Change of velocities in freeze out of QGP to hadronic matter described by massless cut Jiittner distribution. The final 
velocity parameters (full lines) of the cut Jiittner distribition are plotted versus the initial flow velocity of QGP measured in the Rest 
Frame of the Front (RFF) for cases: (a) n =1.2 fm~ 3 , T = 60 MeV, A s = 225 MeV, (b) n =0.1 fm~ 3 , T = 60 MeV, A s = 80 MeV 
(c) no =1.2 fm — 3 , To = 60 MeV, Ag = MeV. Observe that for small initial flow velocities the center of the cut Jiittner distribution 
moves backwards, although all the particles which are allowed to freeze out move forward. Thus, the post FO baryon flow velocities 
(dotted lines) are positive. Note the large acceleration caused by the released latent heat in cases (a) and (b). Ag = B 1 / 4 , where B is 
the bag constant. 



In the general case the solution can be obtained numerically. In the m — limit the solution is simpler and gives 
an interesting insight into the problem. The continuity equation leads to the equation 



Q± 2 (v + l) 3 +v - 1 = 0, 



(18) 



where Q 1 = Q 1 1 (/i,T) — n(/i. T)/(4no7oUo) , which leads to a 3rd order equation and can be solved for v 
alnalytically. The energy equation, [T° M <i(T Al ] = 0, leads to the same equation but with another coefficient 
= Q% (ji,T) — 3Tn(/i, T)/(4eo7ofo) 5 thus, these two equations can have one and the same solution for v, 
only if the two coefficients, Q± and Q2, are equal, which results in 

T = if° 
3 n ' 

and the solutions of both 3rd order equations yield the same expression: 



v = v 3rd (n) = Q 



= O 2 / 3 



l+\/l+Q 2 /27 



1/3 



i-\/i+Q727 



1/3 



- 1. 



(19) 



Then, dividing the equation [T ^d<T M ] = 0, by the equation [T xfJ, da )J ] = 0, yields another 3rd order equation for v: 

R v 3 + 3v 2 + 3(2 - R )v + 3 - 2R a = , 

where R$ = erji>o/po- This equation can be solved analytically and yields one physical root, v — 2 — 3/Rq (and two 
unphysical ones v = —1). Inserting v then into eq. (|l8|), we can obtain the resulting chemical potential, (1, also. 
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Fig. 3. The final baryon density, n, as a function of the pre FO baryon density. The baryon density decreases in the freeze out 
process for cases v = 0.5, T = 50 MeV, (a) A B = 80 MeV (full line), (b) A s = 120 MeV (dashed-dotted line) and (c) A B = 160 MeV 
(dashed line). 

The possibility of this simple analytic solution is a consequence of the fact that in the m = limit the cut of the 
Juttner distribution is made along central cones in the RFG, which then divide the energy and the baryon charge 
exactly in the same proportions. 

As an illustration we studied the freeze out of quark-gluon plasma (QGP) to cut Juttner gas, in-,the massless 
limit. The pre FO side QGP is described by the most simple bag-model EoS (eqs. (5.28-33) in ref.Ea), thus local 
equilibrium is assumed and all pre FO parameters are assumed to be known including the baryon, energy-momentum 
and entropy currents. 

On the post FO side these currents were evaluated earlier in this section, and the equations arising from the 
conservation laws, (|^) , were solved as presented above. Figure || indicates the change of flow velocity during freeze 
out. Physical solution exists only for positive initial velocities, vq > 0. The velocity parameter (!) of the post FO cut 
Juttner distribution varies from —1 to +1, but the post FO Eckart flow velocity is of course always positive in RFF. 
Thus, the post FO baryon current is also positive in RFF (this is obvious since we do not allow any particle to cross 
the front backwards), and consequently, the pre FO current and vo should also be positive because of the continuity 
equation. For small initial velocities, vq — ► 0, the post FO velocities approach zero also, but for moderate velocities, 
deduced recently from experiments, v = 0.3 — 0.7, the difference between the post and pre FO flow velocities may be 
essential. 

In order to show the effect of these modifications compared to the original Cooper-Frye treatment (where the 
increase of the flow velocity is ignored) we can consider case (a) in Figure 2. The cut Juttner distribution always 
leads to an exponential p t spectrum, but according the new modified treatment starting from vq = 0.2c the post FO 
flow velocity increases to vji ow — 0.4c, while the post FO parameter velocity (which determines the pt spectrum) 
increases to v — 0.6c. This corresponds to an increase of the slope parameter, T s i ope , by 60%! This is due to the 
large latent heat arising from the large value of the bag constant taken in case (a). In case (b) the the same effect is 
present but it is weaker. This change of the flow velocity is a basic feature of the correct freeze-out treatment, and 
it is a consequence of the conservation laws and not of the positivity requirement of p^da^ in spacelike FO. Thus, 
the flow velocity change occurs both in spacelike and in timelike freeze-out. This effect can cause for example the 
conversion of latent heat to collective kinetic energy and not to heat if the freeze-out coincides with an exothcrm 
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phase transitional. 

Figure || shows that the baryon density, (^5[), decreases in the freeze out process. This is connected to the fact 
that the post FO flow velocities are above the pre FO ones, as shown in Figure |[ 

We should mention that the post FO temperature parameter of the cut Juttner distribution becomes rather high, 
about an order of magnitude higher than the pre FO temperature. However, we have to recall that the term 
temperature is not applicable for a non-equilibrium distribution, therefore this result has no physical significance, it 
just illustrates the parametrization of the distribution of the assumed cut Juttner sh arae ,.,— ■ 

Finally we have to check the entropy condition for these solutions. As we knowESE3 : E3 QGP can freeze out to 
hadronic matter with entropy production only if the QGP is supercooled or considerably supercooled. This remains 
valid for the cut Juttner assumption as post FO distribution also. With most parametrizations only low temperature 
QGP is able to freeze out. For the cut Juttner gas we cannot speak of a critical temperature, because this gas is not 
in equilibrium and consequently cannot be in phase equilibrium either. Still this distribution can be attributed an 
entropy current by its kinetic definition, and the entropy condition can be checked (Fig. ||). 
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Fig. 4. The ratio of post FO and pre FO entropy currents transverse to the freeze out front. Freeze out can be physically realized 
if R > 1. The entropy condition is tested for three cases: (a) no = 0.1 fm — 3 vq = 0.5, = 80 MeV (full line), (b) no = 0.5 fm — 3 
v = 0.5, A s = 80 MeV (dashed line) and (c) n = 1.2 fm~ 3 v = 0.5, A B = 225 MeV (dotted line). 

In reality the entropy condition is not so stringent as Figure ^ indicates. In this illustrative study the post FO 
EoS had relatively few degrees of freedom to accomodate the high entropy content of QGP. By including many post 
FO mesons and other hadronic degrees of freedom in our post FO EoS, the entropy condition can be satisfied in a 
much wider range of parameters. 




5. Freeze out distribution from kinetic theory 



We have seen that taking the cut Juttner distribution as an ansatz for the post FO distribution, we can solve the freeze 
out problem formally. Although we can satisfy all requirements, the obtained parameter values make it questionable 
whether the cut Juttner ansatz is an adequate assumption. The shape of the distribution with the sharp cuts is also 
a rather unphysical feature of the distribution. 
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To obtain a more realistic, and physically better applicable FO distributions, we should evaluate the distribution 
in more physical, microscopic non-equilibrium models. Kinetic theory is a straightforward candidate for this task—-. 

A first very simplified attempt to solve the freeze out problem dynamically in one dimensional kinetic modelcll 
returned the cut Juttner distribution also, but only in highly unrealistic situations: only when the model yielded 
incomplete freeze out. Thus, further work is needed to find physically realistic post freeze out distributions in kinetic 
models or in other dynamical microscopic models. 

6. Conclusions 

The importance of taking into account conservation laws in the description of the freeze out processes pointed 
out. For freeze out across hypersurfaces with space-like normals the approach suggested by Bugaevj£3 assuming 
cut Juttner distribution as post freeze out distribution is worked out, and the freeze out problem was solved as an 
example for QGP freezing out into a cut Juttner gas. This calculation indicates that results including the Cooper- 
Frye freeze out procedure should be reconsidered and new emphasis should be given to the precise evaluation of the 
post freeze out particle distributions. 

The deviation from the earlier Cooper-Fry approach (where changes of flow velocity, density and temperature were 
ignored) is apparent if the pre FO matter has large energy content in form of compressional energy, latent heat, or in 
any other way, which is not present in the post FO, noninteracting matter. As this post FO matter is not necessarily 
in thermal equilibrium, we cannot consider it as a thermal phase with equilibrium thermodynamical parameters. 
Thus, this idealized approach assuming a FO surface is always assuming a discontinuity irrespective of what was the 
phase of the pre FO matter. Nevertheless, this treatment leads to the strongest modifications in cases when a first 
order phase transition with large latent heat is coupled to the freeze out process. 

Here we have considered an idealized transition as a discontinuity across a hypersurface. In as much as the flow 
across the surface is stationary our results are valid irrespective of the surface thickness, because we used only 
conservation laws. On the other hand in heavy ion reactions the flow across the surface can be considered stationary 
only if it is 1-2 fm wide. With purely kinetic freeze out this is not a very realistic assumption!!] On the other hand 
rapid hadronizaiion from supercooled QGP may satisfy the required conditions and the sharp surface approximation 
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